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Abstract —Approximate expressions for the spatial correlation 
of cylindrical and uniform rectangular arrays (URA) are derived 
using measured distributions of angles of departure (AOD) for 
both the azimuth and zenith domains. We examine massive 
multiple-input-multiple-output (MIMO) convergence properties 
of the correlated channels by considering a number of con¬ 
vergence metrics. The per-user matched filter (MF) signal-to- 
interference-plus-noise ratio (SINR) performance and conver¬ 
gence rate, to respective limiting values, of the two antenna 
topologies is also explored. 

I. Introduction 

In order to meet demands for increased system capacity with 
limited spectral resources, there is a broad consensus that this 
can only be achieved via a large increase in system degrees 
of freedom (d.o.f.) Hl-G). Consequently, massive multiple- 
input-multiple-output (MIMO) systems are being investigated 
m-m, where the number of antennas are scaled up by at 
least an order of magnitude relative to current base station 
(BS) deployments. In the case of massive MIMO, fast-fading 
effectively averages out f7)-[|9], simplifying system analy¬ 
sis. Practical implementation of large numbers of antennas 
required for massive MIMO, typically in confined antenna 
array dimensions, however, results in reduced inter-element 
spacing adversely impacting system performance by way of 
spatial correlation tna. Spatial correlation models for wireless 
systems are thus essential for accurate theoretical performance 
analysis and guarantees. 

An approximate spatial correlation model is presented in 
E2 for clustered MIMO channels, deriving closed-form (CF) 
expressions for a uniform linear array (ULA) and a uniform 
circular array (UCA). This model however does not con¬ 
sider the zenith domain, necessary for accurate performance 
analysis. A correlation matrix is derived in EH which con¬ 
siders both the azimuth and zenith domains, from which it 
is shown that the correlation matrix can be written as the 
Kronecker product of each domains correlation matrix. We 
propose to extend the works of El to massive MIMO antenna 
topologies: namely for uniform rectangular array (URA) and 
cylindrical antenna arrays, while incorporating Third Genera¬ 
tion Partnership Project (3GPP) three-dimensional (3D) urban 
macro cell (UMa) environment measured distributions El- 
With the likelihood of smaller cell sizes and increased line-of- 
sight (LOS) propagation for next generation wireless systems. 


we make a small angle approximation to derive correlation 
expressions. 

The contributions of this paper are as follows: 

• We derive approximations for the spatial correlation of 
a URA and cylindrical antenna arrays using measured 
probability density functions (PDFs) of the AOD in both 
azimuth and zenith domains. 

• We show, via simulation, the impact of antenna array 
topologies on the convergence rate of the channel to 
massive MIMO channel properties. This is shown via a 
number of channel convergence metrics. 

• We show, via simulation, the impact of antenna array 
topologies on performance and rate of convergence, to 
respective limiting values, of matched filter (MF) per-user 
signal-to-interference-plus-noise ratio (SINR). 

• We show that antenna arrays introduce a correlation struc¬ 
ture which does not vary between topology. As a result, 
massive MIMO convergence metrics do not show any 
sensitivity to antenna topology, but do show sensitivity 
to the presence of correlation - which destroys the onset 
of massive MIMO properties. 

II. System Model 
A. System Description 

We consider a multi-user (MU) massive MIMO downlink 
(DL) system with M transmit antennas jointly serving a 
total of K single-antenna users. We assume time division 
duplex (TDD) operation with uplink (UL) pilots enabling the 
transmitter to estimate the DL channel. The M x K channel 
matrix, H, is given by El 

H = R^Hjid, (1) 

where H;;d is the M x K independent and identically dis¬ 
tributed (i.i.d.) channel matrix with CA r (0.1) entries, account¬ 
ing for small-scale Rayleigh fading, and R t is the M x M 
spatial correlation matrix. We consider a cross-polarized (x- 
pol) antenna configuration, with the spatial correlation matrix 
modeled via El 

R t = Xpoi 0 R, (2) 

where the MxM matrix R is the co-polarized (co-pol) spatial 
correlation matrix, © represents the Hadamard product and 



(a) (b) 

Fig. 1: (a) Cluster Model, (b) Channel Model Geometry. 


Xpoi is the M x M x-pol matrix given by 


Xpol — ljf/2 ® 


1 


1 


( 3 ) 


where l M / 2 is a x matrix of ones, S denotes the cross¬ 
correlation between the two antenna elements in the x-pol 
configuration and 0 represents the Rronecker product. 


B. Channel Model 

We consider a clustered channel model, shown in Figure [Ta] 
where we show the zenith AOD offset. Ad, of the wavefront 
relative to the mean zenith AOD of the cluster, 9. In this 
scenario, the linear antenna array is located on the z axis. 
Similarly, if the antenna array was positioned on the y axis, 
Figure [Ta] would describe the offset azimuth AOD, Ac), relative 
to the mean azimuth AOD of the cluster, (j). The channel model 
geometry is clarified in Figure IThl 

When large numbers of antennas are present, antenna 
topologies should exploit x,y and z dimensions. In Section 
nn we derive spatial correlation approximations for the two 
antenna array topologies: URA and cylindrical, depicted in 
Figures [2a] and [2b] respectively. In Figure [2a] the URA is 
geometrically positioned on the y,z plane, where adjacent 
antennas are separated by d\ wavelengths on the z axis 
and ^2 wavelengths on the y axis. Adjacent antennas of the 
cylindrical array, in Figure[2b] are separated by a displacement 
of d\ wavelengths on the z axis and located at a radius of 
p wavelengths from cylinder center, with respect to the x, y 
plane. 

For both antenna topologies, we express the spatial corre¬ 
lation as the Rronecker product of the azimuth and zenith 
domain correlations IH 

R = R 0 0R e, (4) 

and derive the corresponding correlation coefficients indepen¬ 
dently in Section [III] Thus, we consider the URA correlation 
matrix as the Rronecker product of an .4-element ULA in the 
z dimension (zenith domain) and a /1-element ULA on the 
x, y plane (azimuth domain), where M = A x B. Likewise, 
we consider the cylindrical array correlation as the Rronecker 




Fig. 2: (a) URA. (b) Cylindrical Array. 


product of an ,4-element ULA in the z dimension and a 13- 
element UCA on the x , y plane. 

The correlation coefficient between the ath and a'th antenna 
in the zenith domain can then by given by mu 

R*(a,«') = J e j[e ^ Ae) - e “' (Ae)] PAe(A9)dA9, (5) 

where 0 a (A0) is the zenith domain phase shift of the ath 
antennas AOD with respect to a reference antenna, p/\g(A6) 
is the zenith domain AOD offset PDF relative to the mean 
zenith AOD of the cluster. Likewise, the correlation coefficient 
between the foth and //th antenna element in the azimuth 
domain is 

R = J J 

x pA<t>{^)PAe^S)dA(j)dAe, (6) 


where <J>b(Ac)>, A0) is the azimuth domain phase shift of 
the foth antennas AOD with respect to a reference antenna, 
Pa<p(A4>) is the azimuth AOD offset PDF relative to the 
mean azimuth angle of the cluster. Note, we assume the 
azimuth and zenith AOD offset PDF are independent, i.e., 
PA0 ,a<?(A^, A6>) = pAp{A(/))p Ag (A9). 

We model the azimuth and zenith offset AOD PDFs from 
measured values described by 3GPP ED. Thus, the azimuth 
AOD offset, relative to the mean AOD, is modeled as a 
Wrapped Gaussian PDF, given by El 


PA<t>(A<l>) 



0 


(A,» + 27ti) 2 

2 <j^ 

za Act> 


A4> € [—7T, 7r) 

otherwise 

(7) 


where cta 4 , is the standard deviation (SD) of A <f>. Similarly, 
the zenith AOD offset, relative to the mean AOD, is modeled 
by a Laplacian PDF, given by El 


PAe{A0) 


_ I V^Afl I 

^el'^l A0 £ [—7T, 7r) 

0 otherwise 


( 8 ) 


where a A e is the SD of A 8 and k = 1/(1 — e v^OT/a-Ae) 

























normalizes the PDF. 
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C. Convergence Metrics 


In order to study the effects of antenna topology on the 
convergence of massive MIMO properties, we consider the 
convergence metrics used in (8). We evaluate the convergence 
of W = by examining a number of well known 

properties of W and a deviation matrix E = W Ik, where 
I k is the K x K identity matrix. Letting Ai, A2,..., Xk 
denote the eigenvalues of W, we consider: A range. Mean 
Absolute Deviation (MAD) and Diagonal Dominance, defined 
respectively as 


A range — A max ("VV) A m i n (W), 

K 


MAD(E) = —^ J2 l E 

i=ij=i 

W 

1 vv 


yh 


Diagonal Dominance = 


Si=i 


(9) 

( 10 ) 

( 11 ) 


<i=i Jfr ' 

These metrics will be evaluated via simulation for a number 
of system scenarios in Section HVl 


III. Methodology 

We now derive approximations for the spatial correlation 
of URA and cylindrical array topologies. In each case, for 
simplicity we assume no antenna mechanical downtilt. As per 
©, we express the spatial correlation as the Kronecker product 
of the azimuth and zenith domain correlations, and derive the 
corresponding correlation coefficients independently. 


A. Uniform Rectangular Array (URA) 

We first consider the correlation in the zenith domain, Rg, 
of the URA, depicted in Figure [2a] and geometrically described 
in Section lll-BI The zenith domain phase shift of the departing 
wavefront from the ath antenna element, relative to a reference 
antenna, can be expressed as ca 

0 a (A9) = kd\acos{6 + Ad), (12) 

where k is the wavenumber. The angle in (TTH ) can be expanded 
with a first-order Taylor series, while assuming A 9 ss 0, to 
give 

cos(9 + AO) « cos(9) — A0sin(0). (13) 

Since p A g{A9) is zero outside the range [—7r,7r), the PDFs 
integration can be taken over [— 00 ,+ 00 ]. From © we then 
have 

/ OO 

gifcdAa-a') [co.(fl)-Afl B m(8)] pAfl ( A 0) dA 0 

-CO 

(14) 

pOO 

_ ^jkd\ (a— a) cos(0) / jkd\ (a—a') sin(0) A6 


V2a Ae 

— ~jkdl( a ~ a ') COS(0) TT 
- C j 


V2a A9 


V 2 A e I 
I a Ad | 


(16) 


where T’ u denotes the Fourier transform evaluated at u> = 
kd\(a — a') sin(0). Solving the Fourier transform in (IT6l) . we 
have the approximation of the correlation coefficient for the 
zenith domain of a URA as 


R 


0(a,a') 


^jkdi (a—a') cos(0) 

1 + [kd\ (a — a') sin(0)] 2 


(17) 


We now consider the correlation in the azimuth domain. 
The azimuth domain phase shift of the departing wavefront 
from the 6th antenna element, relative to a reference antenna, 
can be expressed as mu 

$b(A<^, A 9) = kd 2 bcos{(f) + Af) sin(0 + A 9). (18) 


We can then use a first-order Taylor series expansion, while 
assuming A cj> sa 0, to express (ITSb as 


cos(^>+A</>) sin(0+A0) sa sin(#+A0) [cos(^) — A^sin(^)] . 

(19) 

Since both p Ac f,(A(f>) and p A g(A9) are zero outside the range 
[— it, it), the integration can be taken over [— oo, +oo]. From 
©, the correlation coefficient between two antenna elements 
6 and 6' is then given by 


R, 


<t>(b,b’) 


fk jkd 2 (b—b') sin(0+A0) cos(0) 


/ oo 

e _ ife d 2 (b-b') sin(f?+A0) sin(0)A^^ A^)dA(j) 

-oo 

xp Ag (A0)dA9. (20) 

Considering the integral with respect to A <j>, embedded in ( 1201 ). 
we have 

pOO 

/ e -jkd 2 (b-b')sin(e+A6) S in(<l>)A<f> p A0 (A0)dA</) 

—jkd 2 (b—b f ) sin(0+A0) sin(</>)A</> 


£ 




_qq CS A< j,y/‘ht 

“1 (A^,+ 2 TTi) 2 

e 


_ (A^.+2»i) 2 

e 2 °U dAf (21) 


bV^tt 


( 22 ) 


VJ=-oo 

where the Fourier transform in (l22l > is evaluated at ui = 
kd 2 (b — b') sin(0 + A 9) sin(</>). Evaluating the Fourier trans¬ 
form, one obtains 

jkd, 2 (b—b r ) sin(9) cos((j>) ^ 


R 


< p(b,b') 
oo 


V2a A g 

e 2 (i-~, 1 2 »i) [ J Aj.fcrf 2 (b-b')sin(6) sin(0)] 2 


—OO 
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i ——oo 


















X T u 


j e 2 (i-j 2 ,ri) [^^kd 2 (b-b') cos(6») sin(0)] 2 (A0) 2 | ^3) 


where steps from (122} to ([23} are given in the Appendix and 
the Fourier transform in (123} is evaluated at 


tu = 


for Ad > 0 
e — 

Ci + for A6 < 0 


(24) 


where 


Ci = jkd 2 {b — b') cos (6) 

x cos(</>)--- cr\ikd 2 (b — b') sin 2 (</>) sin(0) 

j + 2iri v 


(25) 


Evaluating the Fourier transform in ( |23} . we have the approx¬ 
imation of the correlation coefficient for the azimuth domain 
of a URA as 


R-A0(f>,b') ~ nr & 
V 4(7 Ad 


^_ ^jkd 2 (b—b') sin($) cos (<f>) 

!CA6) 

x ^ e 2 (i-~j 2 ,i) [er^kd 2 (b-b') sin(6) sin(c/>)] 2 
i=—oo 

— UJ 2 

x e coa W aln w] (26) 


with ui given in ( 1241) . The spatial correlation approximation 
for the URA is then as per ©, with Re and R,, given in ([17} 
and (126} respectively. 


B. Cylindrical Array 

Similar to the analysis for URA, the spatial correlation of 
a cylindrical array, shown in Figure [2b] can be broken down 
into the Kronecker product of the azimuth and zenith domains, 
described in Section lTl-BI The zenith domain phase shift of the 
ath antenna element, 0„(A6), and the correlation coefficient 
between antennas a and a 1 , Re(a, a'), are identical to the URA 
case and are thus given in (112} and (ITT} respectively. In the 
azimuth domain, the phase shift of the 6th antenna element, 
relative to a reference antenna is given by 05} 

<f> b (A</>, Ad) = fcp6cos((</> — (j>b) + A</>) sin(0 + A 9), (27) 


where <f> — (f>b is the angle between the incident ray projected 
onto the x, y plane and the 6th antenna, relative to the the 
circle center. Note that the azimuth domain phase shift in (127} 
is analogous to (IT8} . with fi = d-j and 6 — <j>b = &■ Therefore, 
the spatial correlation approximation for the cylindrical array 
is given in ©, where 

K 


•R ~ h ~.lkp(b-b )sin(6))cos(0-0 6 ) 

iV 0(6> b ) fK C 
V 4(7 as 

e 2(l4j2,i) [°A<!>kp{b-b') sin(e) sin(0-0 i) )] 2 


E 


p j a^^kp^b — b' ) cos(0) sin(0 — <j: 


x e 1 ~^ 27 


(28) 


where uj is given in (l24} . with 

Ci = jkp(b — 6') cos(6>) x 


cos {(j) - (f) b ) - . xkp{b - b') sin 2 (</> - (j> b ) sin(6») , 

7 + l TTl v 

(29) 

and Re is given in ([17} . 

IV. Numerical Results 

We consider the convergence metrics, described in Section 
III-C1 in order to determine how many antennas are required for 
observable massive MIMO properties in a spatially correlated 
environment. System parameters are presented in Table Q] 


Parameter 

Value 

Frequency (GHz) 

2.6 

X-pol parameter, \J6 

0.1 

Azimuth AOD offset PDF, pAoiAd)) 

Wrapped Gaussian 

Zenith AOD offset PDF, p A e(A0) 

Faplacian 

AOD cluster mean, {4>,9} (log 10 ([“])) 

0.7 

AOD offset SD, {o A4> ,o Ae } (log 10 ([°])) 

-0.3 


TABFE I: System Parameters 


A. Convergence Properties 


M=50,URA 
M=50,URA Approx 
M=50, Cylindrical 
M=50,Cylindrical Approx 
M=100,URA 
M=100,URA Approx 
M=100,Cylindrical 
M=100,Cylindrical Approx 
M=200,URA 
M=200,URA Approx 
M=200,Cylindrical 
M=200,Cylindrical Approx 

“0 5 10 15 20 25 30 

X range 

Fig. 3: A range CDF for f- = ct = 10 

In Figure [3] we show the cumulative distribution function 
(CDF) of the A range, given in ([9}, for varying M. The 
eigenvalues of W are generated using H, given in ©, where 
the simulated CDFs are generated from instantaneous array 
factors, while the approximations are generated from ©. It 
is seen that as we increase the number of transmit antennas, 
M, the median value of the A range CDF increases, rather 
than converging to an equal eigenvalued channel, which is 
observed in El and shown in 0 for an i.i.d. channel. This is 
due to an increase in the dominant eigenvalue of W, resulting 
from such a narrow angle spread. From Figure [3] we see that 
for all values of M, our derived expressions approximate the 
spatial correlation well. Thus, we use the approximations for 
all results following. 














































Fig. 4: MAD(E) vs K, with ^ = a = 10 and 25 


In Figure [4] we plot the average of MAD(E), given in ( ITOb . 
versus K for a = 10 and 25. We observe very different 
behaviour between the correlated and i.i.d. scenarios. In the 
i.i.d. case, MAD(E) converges slowly to zero for increasing 
K, with a quicker convergence rate for larger a. Note that the 
average MAD(E) value is plotted in Figure [4] and the only way 
for this to converge to zero is for all W matrices to be close 
to I k- Hence, only in the i.i.d. case does each W become 
close to I k as K increases. 



Fig. 5: Diagonal Dominance vs M, with K fixed and I\ = ^ 

Figure [5] shows the diagonal dominance of W (size K x K) 
given in £□>, as a function of the number of transmit antennas, 
M. We see that the correlated scenario has similar trends 
as the i.i.d. case. As with the results in j 8 ) for K = —, 
the diagonal dominance decays. This follows as the number 
of off-diagonals, which grow at a rate « K 2 , increase the 
denominator of (fill at a faster rate than the numerator. In 
the correlated scenario, the diagonal dominance approaches 
zero much faster than the i.i.d. case. Here, while the diagonal 
elements converge to a mean of 1 , the large number of 


off-diagonals converge to a non-zero mean, dominating the 
diagonal elements. On the other hand, for fixed I\, the sum 
of the diagonal elements increase by a greater proportion than 
the off-diagonal elements as M is increased, resulting in a 
more diagonally dominant i.i.d. W. In the correlated case 
the diagonal dominance converges quickly to small non-zero 
value. 

Considering Figures [3] 0 ] and 0 ] we conclude that while 
for i.i.d. channels, the massive MIMO metrics converge for 
a very large numbers of antennas, desirable massive MIMO 
properties are degraded with spatial correlation present (with 
small angle spreads). 

B. Convergence Properties of MF Precoder 

We now explore the impact of massive MIMO antenna 
topologies on MF SINR performance and convergence to 
limiting values. The MF SINR of the ith user is given by 



where p& is the DL transmit signal to noise ratio (SNR), h, 
denotes the ?th column of H, and 7 = tr(H T H*)/Ff is the 
power normalization factor. 



Fig. 6 : MF SINR for, ^ = a = 10, p d = 10 dB 


In Figure [ 6 ] we plot the expected value of ( 130k It can be 
observed that there is a huge reduction in E[SINRj] perfor¬ 
mance by introducing correlation of which antenna topology 
has almost no effect on. The MF SINR for an i.i.d. scenario 
rapidly converges to its limiting value, whereas in a correlated 
case, E[SINRj] converges to limiting values for M > 200. 
We conclude that linear precoders are highly sub-optimal for 
massive MIMO systems in spatially correlated environments 
with small angle spreads. 

V. Conclusion 

In this paper we develop approximations for the spatial 
correlation of URA and cylindrical antenna arrays, from a 3D 
channel model. Using derived approximations, we show that 























































large spatial correlation, due to small angle spread, destroys 
the convergence of massive MIMO properties with increasing 
number of antennas. Furthermore, the impact of massive 
MIMO antenna topology considered is shown to be negligible. 
MF SINR performance and convergence rate is explored for 
the two antenna topologies, showing the detrimental effect of 
correlation, with respect to an i.i.d. channel. It should be noted 
that if we increased inter-element spacings of the antenna 
arrays relative to the angle spread, we would see a decrease in 
correlation and the convergence metrics would approach the 
i.i.d. case. We leave this to future work. 


VI. Future Work 

In the future, we aim at deriving correlation matrices 
for various antenna topologies which do not rely on small- 
angle distributions. Also, we aim to incorporate the effects of 
mutual coupling, due to the proximity of antennas as electrical 
components, important for the analysis of large antenna arrays. 
Furthermore, an investigation into how the various antenna 
spacing parameters influence the correlation for the antenna 
topologies will give insight into the structure of the correlation 
matrices and how to carefully design an antenna array. 


Appendix 


From ( 1221 . we have 


R, 


4>{b,b') 


f 


jkd 2 (b—b f ) sin(0+A0) cos ((f)) 


e^ 1 -^ 


:[a^kd 2 (b-b')sin(e+Ae)sin( C l>)}' 2 pAg ^ A ^ dAd ( 3 ^ 


~jkd 2 (b— b / )[sin(0)+A0 cos(0)] cos(0) 


E e’u-A") 


X [cr A ^fcd 2 (b-b') [sin(6»)+A0 cos(0)] sin(0)] 


_ ^jkd 2 {b—b') sin(0) cos(0) 


E 


e 2 (l-j 2 „i) 


X \a^kd 2 (b— b') sin(0) sin(0)] 2 / e jkd 2 (b—b') cos(S) 

J —oo 

X [cos(0)~ o- A</ ,fcd 2 (b-b') sin 2 (0) sin(0)]A0 x 

e 2 (i~j 2 -Ki) K *kd 2 (b-b') cos(0) sin(0)] 2 (■ Ae )' 2 pAg (A6») dA9 (32) 

OO 

— c ikd 2 (b-b') sin(6>) cos (<f>) K e 2(i-/ 2 7ri) 

V^cTAe ■ 

v i= — oo 

X [&A(t>kd 2 (b— b') sin(0) sin(0)j / ^jkd 2 (b—b') cos(Q) 

J — ( 


x[cos(0)—j T U FT CT A</ ,fcd2(b-b')sin 2 (0)sin(0)]A6( x g-^jl Ae l 

x e 2 (i-jLi) [^A^fcd2(b-b')cos(0)sin(0)] 2 (A9) 2 ^^g ^ 33 ^ 

where a Taylor series expansion is used to obtain (|32T >. 
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